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Methods for distilling maximally entangled tripartite 
(GHZ) states from arbitrary entangled tripartite pure states 
are described. These techniques work for virtually any input 
state. Each technique has two stages which we call primary 
and secondary distillation. Primary distillation produces a 
GHZ state with some probability, so that when applied to an 
ensemble of systems, a certain percentage is discarded. Sec- 
ondary distillation produces further GHZs from the discarded 
systems. These protocols are developed with the help of an 
approach to quantum information theory based on absolutely 
selective information, which has other potential applications. 

In the rapidly developing field of quantum informa- 
tion, it is possible to identify two main lines of investi- 
gation. On the one hand, it addresses basic questions on 
the fundamental nature of information, how it is embod- 
ied in quantum systems, how it can be quantified, and 
the extent to which physical properties can be reduced to 
informational ones JjJ . On the other hand are specific op- 
erational issues; for example, how quantum information 
can be manipulated for applications such as quantum 
computation and teleportation 

In this Letter we try to bring together these two 
strands by proposing a new approach to the analysis 
of quantum information at a fundamental level, which 
leads directly to an operational technique for distill- 
ing maximally entangled tripartite (GHZ) states ||, us- 
ing local operations and classical communication. The 
three qubits of the system are assumed to be physically 
separated and held by Alice, Bob, and Cara, respec- 
tively. (Throughout this Letter we use the term "maxi- 
mally entangled" to refer to N-partite states that are N- 
orthogonal: i.e., if the subsystems are two-dimensional 
such a state would be 072(10000 . . .) - |1111 . . .)).) 

Central to our approach is the notion of absolutely se- 
lective information, which has a straightforward interpre- 
tation in terms of classical information but can be seen as 
a basic distinguishing feature between quantum systems 
and their classical counterparts. We apply our approach 
to the specific problem of distilling GHZ states from ar- 
bitrary entangled tripartite pure states. We show that it 
is possible to distill, with a certain probability, a GHZ 
state from virtually any entangled tripartite pure state 
while retaining all three subsystems of the input state. 
As far as we know this is the first protocol of this type 
to be suggested. Our initial yield of GHZ states is then 
supplemented by an additional yield which involves sac- 



rificing some subsystems. In this Letter we outline our 
approach and summarize our results. A more detailed 
exposition of the underlying analysis will be presented 
elsewhere H. 

The distinction between "selective" and "structural" 
information was addressed by Mackay || in the early 
days of classical information theory. Whereas structural 
information measures are based on an analysis of the 
form of possible events, selective information refers to 
new information gained from the occurrence of a specific 
event. For example, a signal might transmit a bit as 
one of two different waveforms; the selective information 
would be one bit, while the structural information, suf- 
ficient to describe all possible waveform measurements, 
would be considerably more. Absolutely selective infor- 
mation signifies data that are irreducibly unpredictable, 
and hence genuinely new, in the sense that their unpre- 
dictability cannot be explained by the observer's igno- 
rance. This type of information can arise only in a the- 
ory that is fundamentally stochastic; hence it is com- 
monplace in quantum physics, but absent from classical 
physics. For a pure state, the minimum local absolutely 
selective information (the minimum information gener- 
ated by measuring one of the subsystems with a free 
choice of measurement basis) is exactly the same as the 
local entropy. 

When considered as a quantitative measure, selective 
information is closely related to fundamental measures 
in quantum information theory. For example, the stan- 
dard measure of entanglement for bipartite pure states 
Q is numerically equal to the minimum local absolutely 
selective information. In a similar way, minimizing the 
absolutely selective information can be used to develop 
measures of nonorthogonality for quantum states . In 
this Letter we show that absolutely selective information 
can be manipulated by an appropriate measurement pro- 
cedure, and apply this to an operational problem. 

The problem we address is to transform a state \1p123} 



^123) = a|000) + 6|001) + c|010) + d|011) 

+e|100) + /I101) + $|110) + &|111) 



(1) 



into the state |V>ghz) = y/T/2(\000) - I 111 )). witn somc 
probability. Let us first consider the minimum absolutely 
selective information Ai associated with each of the three 
qubits in ^123): 
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Ai = - (pilog 2 Pi + (1 -pi)log 2 (1 % = 1,2,3, 

(2) 

where pi and 1—pi are the eigenvalues of the reduced den- 
sity operator pi that describes system i when the other 
two subsystems are traced out, 



Pi = Tr jk (\^ijk){ipijk\) 



(3) 



(Without loss of generality we adopt the convention that 
Vl > 1/2.) 

It can be shown [g that the maximal value of Ai 
for any tripartite pure state occurs uniquely for the GHZ 
state (or any local unitary transform of it), for which 
each of the pi's is equal to 1/2 and Yli=i Ai = 3. Our 
distillation procedure makes use of this fact by decreas- 
ing each of the p^'s in turn (with some probability), ap- 
plying the procedure repeatedly, until all of the p^s are 
within some tolerance of 1/2, at which point a GHZ state 
will necessarily have been distilled (up to a local unitary 
transformation) . 

The technique for decreasing the pj's is similar to the 
"Procrustes" method of ||. We perform a positive op- 
erator valued (POV) measurement consecutively on each 
subsystem. To see how this works, let us carry out the 
procedure on subsystem 1. If we consider subsystems 2 
and 3 as a single composite system, we can write the tri- 
partite state as a Schmidt decomposition with respect to 
subsystem 1 and the composite 2-3 system: 



123/ 



pri+)il^ + >23 + v / i 



Pi 



/23- 



(4) 



The minimum absolutely selective information for sub- 
system 1 is then given by 



Ai = - (pilog 2 pi + (1 -pi)log 2 (1-pi)) 



(5) 



By carrying out an appropriate POV measurement on 
this subsystem, we can with some probability bring Ai 
to its maximal value of 1. We introduce an ancilla qubit 
"a", which interacts unitarily with subsystem 1: 



l+)i|0) Q 

|-)ilQ) 



«i+)iio> a + v / r^i+> 1 ii> Q , 

l-)ilo> a , 

Vl-a2| +)i | ) a -a|+) 1 |l) a , 
Hill).- 



(6) 



We then measure the state of the ancilla. If we set 
a = y/ (1 — pxf/px and the starting state of the ancilla 
to be |0) , we will with probability 2(1— pi) find the 
ancilla in state |0) , which projects the system into the 
state yi72(|+) 1 |0 + ) 23 + h)il0~) 23 ), for which A 1 = 1. 
With probability 2pi — 1 we will measure state |l) a , in 
which case the procedure fails. 

In the Procrustes technique, this one step plus a local 
unitary transformation suffices to distill EPR pairs from 



arbitrary entangled bipartite states ||. In the tripar- 
tite case, we then repeat the procedure on subsystems 2 
and 3, which projects the system into states for which 
p 2 = 1/2, A 2 — 1 and p% = 1/2, A 3 = 1, respectively. 
However, if we simply carry out a single POV measure- 
ment of this type on each of the three subsystems in 
turn, the resulting tripartite state will not in general be a 
GHZ state. Each step of the process is nonunitary, since 
the tripartite system can be discarded at each stage if 
the wrong result for the POV measurement is obtained. 
Whilst the A4 are conserved by unitary operations on 
the local subsystems, they are not conserved in general 
for nonunitary operations. Hence, when we carry out a 
POV measurement on bit i to project the system into a 
state for which pi = 1/2 and Ai — 1, this will disrupt the 
values of p and A for the other two qubits. 

Nevertheless, it transpires that for most tripartite 
states, repeated application of this type of POV mea- 
surement will steadily move the input state towards a 
GHZ state until it gets arbitrarily close to it. (Excep- 
tions will be identified later). There are a number of 
plausible ways to measure "closeness" to a GHZ state. 
Three such measures are 



D p =J2 Pi -3/2, 
i=l 

3 

D s = 3-Y,Ai, 

i=l 

3 

D 2 = 3/4 -$>i(l -Pi). 



(7a) 
(7b) 
(7c) 



We introduce this last quantity because it is more 
tractable analytically than D p and D$, being a simple 
function of the coefficients of |V>i23)- 

Numerical analysis of this process on randomly chosen 
starting states shows that for a large fraction of these 
states, Dp approaches zero to an accuracy of 10~ 3 after 
just two complete iterations (i.e. two POV measurements 
performed on each of the three subsystems), while virtu- 
ally all do so within four iterations. Interestingly, we 
find that in every case examined (aside from the excep- 
tions given below), D p decreases monotonically toward 
zero with each step of the procedure, whereas D$ and 
.D 2 can fluctuate, though of course their general trend is 
downward. 

The results presented so far are supported only by 
numerical analysis. However, there is a closely related 
procedure for which we have derived an analytical proof 
of efficacy for virtually any input state. In this second 
method, instead of reducing each probability pi to 1/2 
in turn using POV measurements, the probabilities are 
reduced by a small amount e, so that with each step 
the state changes infinitesimally in the limit e — * 0. 
The proof follows fairly straightforwardly by deriving the 
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changes in the state coefficients from the procedure, then 
using these to derive the change in £>2- By changing to 
the Schmidt basis for all three bits, one can, with some ef- 
fort, show that pi(l—pi) can never decrease for i = 1,2, 3, 
and will only remain unchanged for certain very special 
initial states detailed below. This monotonicity implies 
monotonicity for D2, D p , and D$, so all of these quan- 
tities diminish steadily as the state approaches a GHZ 
state. This infinitesimal method would be quite chal- 
lenging experimentally, but it is analytically interesting 
due to its relative tractability. 

The protocols described so far correspond to what we 
call "primary" distillation. They will give a specific yield 
(i.e., surviving percentage) of GHZ states if a collection of 
systems is supplied in a given input state. This yield can 
be straightforwardly calculated for the large-step proce- 
dure; after each POV measurement on the ith subsystem 
a proportion 2(1 — Pi) of the systems are retained. The 
yield of GHZ states for the primary distillation process 
averaged approximately 9.2% for the evenly-distributed 
sample of input states we analyzed, but this will clearly 
depend strongly on the initial distribution. 

Average yields for the infinitesimal procedure were 
9.7%. The chance of the procedure failing on any given 
step is quite small, but over many steps the number of 
failures mounts. The difference between the infinitesimal 
and big-step procedure is interesting when contrasted 
with the bipartite Procrustes technique. In the bipar- 
tite case, there is no advantage to using small steps over 
a single large step; the yields are the same in both cases. 
Clearly in the more elaborate tripartite procedure there 
is a difference. 

This yield can be greatly enhanced by a process of 
secondary distillation, which makes use of those systems 
discarded during primary distillation. When we carry 
out the initial POV measurement on subsystem 1 for 
the input state IV'123) given by eq. (0), with probability 
2pi — 1 we will fail to obtain the desired result. However, 
this failure will leave the discarded system in the state 
|+) 1 |</> + ) 23 , where |</> + ) 23 is in general an entangled bi- 
partite state of subsystems 2 and 3. Similarly, failures at 
later steps of the primary distillation process can yield 
entangled bipartite states of subsystems 1 and 2 and of 
subsystems 1 and 3. Thus, when the primary distillation 
procedure has been completed on a collection of systems 
in a given input state, we will have an additional residue 
of entangled bipartite states of subsystems 1 and 2, 1 and 
3, and 2 and 3. These entangled pairs can be distilled to 
EPR pairs by standard techniques (6), and the resulting 
EPR pairs can be used to prepare further GHZ triplets. 

This is quite similar to the method of Q , where GHZ 
states were produced by first distilling EPR pairs be- 
tween Alice, Bob, and Cara, and then using two pairs 
to produce a GHZ triplet. (For example, if Alice shares 
one EPR pair with Bob and another with Cara, she can 
distribute a GHZ state by preparing it locally and then 



teleporting the states of two of the subsystems to Bob 
and Cara with the help of the two EPR pairs.) If, when 
primary distillation is completed, we produce A23 EPR 
pairs of subsystems 2 and 3, A31 EPR pairs of subsystems 
3 and 1, and N12 EPR pairs of subsystems 1 and 2 from 
the discarded systems, we will be able to distill a further 
(A23+A31+A12V2 GHZ triplets (in the case where none 
of the As is greater than the sum of the other two), or 
(Njk + N ki ) GHZ triplets (if > N jk + N ki ). Numeri- 
cal analysis indicates that the average yield for secondary 
distillation of GHZ states, for the random sample consid- 
ered, is approximately 27.5% giving a total yield of about 
36.7%. The infinitesimal technique does even better, giv- 
ing a secondary yield of 29.3% for a total yield of 39%. 

Since the bulk of this yield comes from the production 
of EPR pairs, one might reasonably ask how these meth- 
ods compare to simply producing EPR pairs (with no 
primary distillation), and then using these pairs to pro- 
duce GHZ triplets directly EPR pairs are produced 
by measuring one of the subsystems in such a way as to 
maximize the pair-wise entanglement between the other 
two bits, and then distilling perfect EPR pairs from the 
resulting states. For the same random sample of states, 
this technique produces an average yield of 31.5%, lower 
than either of the other two techniques, and not much 
higher than the secondary yield alone. This does not, of 
course, prove that it is worse for every initial state. How- 
ever, the closer the initial state is to a GHZ (using any 
of our distance measures (J7a|--|7c|) ) , the better these distil- 
lation procedures perform, while producing GHZs from 
EPR pairs has a maximum yield of 50%. (It may be pos- 
sible to increase this yield asymptotically; but if, so the 
secondary yield in our protocols will also be increased, so 
we do not expect this to change our conclusions substan- 
tially.) 

For some special cases these protocols will not work 
as described. If the original input state is not three- 
party entangled, the protocol will fail completely; that 
is, if the original state can be written as Ix^lOjfci no 
three-party entanglement will be distillable by either pri- 
mary or secondary distillation. There is another set 
of states for which primary distillation fails, but which 
can still produce GHZ states by secondary distillation. 
This set consists of tripartite input states with just 
three components, where each component is biorthog- 
onal (but not triorthogonal) to the other two, and local 
unitary transforms of such states. For example, the state 
\ip tT ) = b|001) + c|010) + ejlOO) is of this type. We call 
such states "triple" states; all have D p > 1/2, though a 
substantial subset has D p = 1/2 exactly. Both forms of 
the primary distillation process take triple states to triple 
states, so that the GHZ state will never be produced. The 
large step procedure causes all triple states to converge 

to the state \ipau) = y/T/2\001) + J (y/Z - 1)/4|010> + 
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J (3 - \/5)/4|100), at which point any further steps will 
simply result in a cyclic shuffling of the component am- 
plitudes. We call this attractor state the "golden mean" 
state because of the appearance of the golden mean in 
the amplitudes. The infinitesimal procedure leaves all 
triple states with D p = 1/2 unchanged. This procedure 
can also cause certain other states with D p > 1/2 to 
converge to triple states rather than GHZ states, though 
most do not. The large step procedure may also take 
some states to triple states. Although triple states do 
not yield any GHZ states by primary distillation, they 
can of course produce them by secondary distillation. 

What is more, it is possible to move off of a triple 
state (with some probability) by performing a POV 
measurement in a basis other than the Schmidt basis. 
That is, instead of using the basis and in the 
transformation (H), one uses a different basis, such as 
± |— ) 1 )/v / 2^ Setting a to a reasonable value (such 
as a — \/l/2) will then take triple states to non-triple 
states. In principle, the state might then re-converge to 
a triple state, but this is quite unlikely. 

We have shown that manipulation of absolutely selec- 
tive information can be used to distill maximally entan- 
gled tripartite states from arbitrary tripartite entangled 
pure states. This method will not work for systems with 
four or more subsystems, since in these pi — 1/2 does 
not uniquely determine the maximally entangled state. 
It may be that related techniques might succeed, how- 
ever, if it is possible to manipulate other locally unitarily 
invariant parameters by local POVs and classical com- 
munication. Even in the 3 qubit case, however, the pro- 
cedure we have described is surely not optimal. The op- 
timal distillation technique is not known, but would al- 
most certainly make use of joint manipulations on many 
copies of the input state. Nor is an asymptotically re- 
versible distillation technique known for GHZ states jlO| . 
It would be interesting to compare the yields of these 
two hypothetical techniques. In the bipartite case they 
are the same, but this need not be so in the tripartite 
case. Indeed, if the reversible GHZ distillation technique 
produced any extra two-party entanglement, one would 
generally expect to be able to produce further GHZ states 
by an irreversible secondary distillation stage. This sug- 
gests that the algorithm giving the optimal yield of GHZs 
will probably not be reversible. It would also be interest- 
ing to compare our yield of GHZ states to some standard 
measure of tripartite entanglement. Lacking such a mea- 
sure, however, the best that can be done is to compare 
different distillation techniques to each other. 

There may be a number of other problems in quantum 
information theory which are amenable to an approach 
focusing on the absolutely selective information content 
of quantum systems. For example, work in progress sug- 
gests than such an approach can be useful in the analysis 
of nonorthogonality. Since selective information is a clas- 



sical concept, this approach also provides a valuable link 
between classical and quantum information. 
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